AD-A268  650 


AIAA-93-0631 

Corrections  to  Fringe  Distortion  due  to 
Flow  Density  Gradients  in  Optical 
Interferometry 

Y.C.Cho,  NASA  Ames  Research  Center, 
Moffett  Field,  CA; 

L. W.Carr,  U.S.Army  ATCOM, 

NASA  Ames  Research  Center, 

Moffett  Field,  CA;  and 

M. S.Chandrasekhara,  Naval  Postgraduate 
School,  Monterey,  CA 


DTIC 

ELECTE  I 
AUG  3  01993  I 


This  document  has  been  approved 

tor  public  release  and  sale;  its  QQ— 

distribution  is  unlimited. _ 

- o  n  9  ■■■' 

31  st  Aerospace  Sciences 
Meeting  &  Exhibit 

January  11-14, 1993  /  Reno,  NV 


93-19880 

liiiimiH  \\pf 


For  permission  to  copy  or  republish,  contact  the  American  Institute  of  Aeronautics  end  Astronautics 
370  L*Enfant  Promenade,  S.W.,  Washington,  D.C.  20024 


THIS  DOCUMENT  IS  BEST 

I.. 

QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


Corrections  to  Fringe  Distortion  due  to  Flow  Density 
Gradients  in  Optical  Interferometry 

Y.  C.  Cho1 

Information  Sciences  Division 

NASA  Ames  Research  Center,  Moffett  Field,  CA  94035-1000 


'3 

3 


L.W.  Carr2 

Aerofiightdynamics  Directorate,  U.S.Army  ATCOM 
and  Fluid  Mechanics  Laboratory  Branch 
NASA  Ames  Research  Center,  Moffett  Field,  CA  94035-1000 

M.S.  Chandrasekhara3 
Navy-NASA  Joint  Institute  of  Aeronautics 
Department  of  Aeronautics  and  Astronautics 
Naval  Postgraduate  School,  Monterey,  CA  93943 

DTIC  qUAUTY  CIEP2CT3D  3 


Abstract 

An  analytical  method  is  formulated  to  account  for 
distortions  of  optical  interferograms  used  for  studies 
of  flow  over  airfoils  experiencing  dynamic  stall.  It  is 
shown  that  such  distortions  are  generated  primarily 
due  to  optical  path  deflections  in  the  test  flow  ,  caused 
by  large  density  gradients.  Such  in-flow  optical  path 
deflections  are  neglected  in  conventional  optical  tech¬ 
niques  for  flow  studies.  The  present  method  employs 
a  ray  analysis  to  determine  these  in-flow  optical  path 
deflections,  and  accurately  predicts  the  interferogram 
distortions 

Nomenclature 

a  coefficient  of  y  in  expansion  of  air  density 
as  a  function  of,  Eq.  (11) 

b  coefficient  of  y2  in  expansion  of  air  density 
as  a  function  of  y,  Eq.(U) 

e  base  of  natural  logarithms 

/  focal  length  of  spherical  mirror,  Fig.  3 

I  light  intensity 

Jo  intensity  of  an  incident  light 
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K  Gladstone- Dale  constant,  Eq.  (1) 

L  wind  tunnel  span  (25.4  cm),  Fig.  2 

£  Thickness  of  glass  window,  Fig.2 

M  Free  stream  Mach  number 

m  image  magnification,  Eq.  (25) 

n  (spatially  varying)  index  of  refraction  of 
air 

p  distance  from  mirror  to  an  object  plane, 
Fig.  3 

q  distance  from  mirror  to  an  image  plane, 
Fig.  3 

s  one-dimensional  spatial  coordinate 

variable  along  a  ray 

w  a  parameter,  Eq.  (15) 

z,  y,  z  spatial  coordinates:  x  -  mean  flow 
direction,  y  -  vertical,  z  -  spanwise 

Vo  y-coordinate  of  an  incident  ray  at  z=0 

yt  y-coordinate,  Eq.  (26)  and  Fig.4 

PL  y-coordinate  of  a  ray  at  exit  plane 

z=L 

tangential  slope  of  a  ray,  Eq.  (17) 
tangential  slope  of  a  ray  at  exit  plane  *=L 
a  angle  of  attack  of  airfoil 

0  y/2Kb,  Eq.  16 

7  Ap/po 

6  virtual  displacement  at  an  object  plane, 

Eq.  (22)  and  Figs.  2  and  3 

60  virtual  displacement  at  s=0 

61  virtual  displacement  at  x=L 

6  image  displacement,  Eq.  (23) 
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So  image  displacement,  focused  on  object 

plane  z=0 

6 1  image  displacement,  focused  on  object 

plane  z=L 

9  angle  of  deflected  ray  measured  from 
z-axis 

p  (spatially  varying)  air  density 

po  ambient  air  density 

p,  air  density  at  surface  of  airfoil 

poo  air  density  in  free  stream 

A  p  poo  -  p, 

1.  Background 

A  real-time  optical  interferometry  technique  has 
recently  been  used  at  NASA  Ames  Research  Center  to 
study  dynamic  stall  over  an  oscillating  airfoil  in  com¬ 
pressible  flows.1,2  During  this  investigation,  a  time 
sequence  of  interferograms  was  documented  and  ana¬ 
lyzed.  Each  interferogram  typically  displays  a  fringe 
pattern  delineating  the  density  field  surrounding  the 
airfoil.  The  separation  of  two  adjacent  fringes  cor¬ 
responds  to  a  phase  increment  of  2x.  The  fringe  in¬ 
tensity  varies  sinusoidally  along  a  line  passing  from 
one  fringe  to  the  next  one,  implying  a  linear  change 
of  the  phase  and  of  the  flow  density  within  the  fringe 
separation.  As  the  flow  velocity  increases,  the  fringe 
frequency  increases  locally  in  regions  near  the  airfoil 
surface,  indicating  the  presence  of  a  large  density  gra¬ 
dient. 

At  high  Mach  numbers,  the  interferogram  suffers 
distortions  locally  where  the  fringe  frequency  is  very 
high.  Evidence  of  such  distortions  is  exhibited  by 
fringes  that  separate  from  or  penetrate  into  the  image 
of  the  upper  surface  of  the  airfoil  at  angles  of  attack  a 
=  5“  ~  15°,  and  flow  velocities  M  =  0.3  ~  0.45,  that 
is,  when  the  airfoil  is  experiencing  very  high  levels 
of  lift  [see  Plate  I],  Another  example  of  the  distor¬ 
tion  is  manifested  near  the  leading  edge  of  the  airfoil 
for  large  flow  velocities.  In  this  case,  the  fringes  dis¬ 
appear  locally  leaving  the  interferogram  with  a  com¬ 
pletely  dark  spot  as  shown  in  Plate  II.  Detailed  under¬ 
standing  and  correction  of  interferogram  distortions  is 
critical  to  accurate  analysis  of  the  compressible  flow 
that  occurs  on  airfoils  at  high  angles  of  attack.  Cor¬ 
rection  of  these  distortions  is  absolutely  necessary  if 
quantitative  analysis  of  the  flow  is  to  be  considered 
in  this  unsteady  flow  environment. 

The  primary  cause  of  the  observed  fringe  distor¬ 
tions  is  evidently  optical  refraction  induced  by  the 
large  density  gradients.  An  optical  refraction  in  a  test 
flow  may  be  regarded  as  composed  of  two  elements: 
Phase  modification  and  path  deflection.  Although  the 
in-flow  deflections  are  effected  as  modifications  of  the 
wave  front  defined  as  a  surface  of  constant  phase,  it 
is  neglected  in  conventional  optica]  methods  applied 
to  flow  studies  such  as  shadowgraphs  and  schlieren 
systems.  These  methods  assume  that  the  deflection 
occurs  only  at  the  interface  of  the  test  flow  and  the 


surrounding  medium.  Even  in  optical  interferometry, 
the  in-flow  deflections  were  previously  neglected  be¬ 
cause  their  effects  were  insignificant  or  appeared  to 
be  insignificant.  We  initially  took  such  neglect  for 
granted,  and  suspected  that  a  potential  cause  of  the 
observed  fringe  distortions  were  the  double  refractions 
of  the  light  at  both  surfaces  of  the  exit  glass  window. 
However,  an  in-depth  investigation,  both  analytical 
and  experimental,  showed  that  such  double  refrac¬ 
tions  would  displace  image  planes  but  could  not  gen¬ 
erate  the  observed  image  distortions.  In  this  paper, 
an  analytical  method  has  been  formulated  to  account 
for  the  interferogram  distortions  including  in-flow  op¬ 
tical  path  deflections  which  are  determined  by  using 
a  ray  analysis. 

2.  Experimental  Setup 

In  the  present  study,  the  airfoil  is  mounted  in  a 
wind  tunnel  between  two  2.54  cm  thick  glass  windows 
which  are  oriented  perpendicular  to  the  spanwise  co¬ 
ordinate  z-axis  and  parallel  to  each  other.  The  in¬ 
ner  surfaces  of  the  two  glass  windows  are  specified  in 
terms  of  two  planes  z  =  0  and  z  —  L  (=25.4  cm).  The 
flow  is  confined  between  these  two  planes,  and  is  as¬ 
sumed  to  be  two-dimensional;  thus  the  field  quantities 
-  the  pressure,  the  density,  and  the  flow  velocity,  are 
independent  of  z.  The  interferograms  were  obtained 
using  a  Point  Diffraction  Interferometry  technique3 
using  a  Nd:YAG  laser  operating  at  the  optical  wave¬ 
length  532  nm  (see  Fig.  1).  The  laser  beam  is  colli¬ 
mated  so  that  all  the  light  rays  passing  through  the 
flow  are  parallel  to  the  z-axis.  The  interferogram  is 
obtained  as  a  sinusoidal  function  of  the  phase  dif¬ 
ference  between  the  reference  beam  and  the  signal 
beam  that  experiences  the  phase  modification  during 
its  passage  through  the  flow.  The  reference  beam  here 
also  passes  through  the  flow  and  the  optical  imaging 
system  but  without  any  phase  modification.  Thus, 
the  point  diffraction  interferometry  is  immune  to  in¬ 
terferogram  modifications  due  to  diffraction  (Fres¬ 
nel  or  Fraunhofer)  encountered  in  usual  holographic 
interferometry4,6. 

3.  Ray  Analysis  for  Optical  Path  Deflections 
in  Test  Flow  with  Large  Density  Gradients 

As  mentioned  earlier,  the  observed  fringe  distor¬ 
tions  are  caused  by  the  optical  refraction  induced  by 
the  large  density  gradient.  However,  the  conventional 
treatment  of  the  refraction  (neglecting  the  ray  deflec¬ 
tion  in  the  test  section)  is  not  adequate  to  account  for 
the  fringe  distortions  observed  in  the  present  exper¬ 
iment;  thus  the  imaging  has  been  more  thoroughly 
investigated  by  using  an  analysis  that  includes  ray 
deflection  in  the  test  section.  In  doing  so,  we  con¬ 
sider  a  simple  flow  in  which  the  air  density  in  the  test 
section  varies  only  in  the  y-direction  (see  Fig.  2  for 
coordinates).  In  reality,  the  air  density  varies  in  the 
x-direction  as  well  as  in  the  y-direction.  The  analy¬ 
sis  can  be  easily  applied  to  this  problem  by  replacing 
the  y-coordinate  variable  with  a  one-dimensional  co¬ 
ordinate  variable  which  is  perpendicular  to  the  t-axis. 
The  new  coordinate  axis  is  in  the  direction  of  Vp,  in  a 
cylindrical  region  within  which  a  ray  under  investiga¬ 
tion  is  confined.  In  the  two-dimensional  flow  as  in  the 
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present  study,  the  cylinder  is  parallel  to  the  z-axis. 

The  density  variation  causes  variation  of  the  re¬ 
fractive  index  in  accordance  with  the  Gladstone- Dale 
equation,6 

n(y)  =  1  +  Kp(y).  (1) 

where  K  is  the  Gladstone- Dale  constant  whose  value 
is  .227  cm3/g  for  the  optical  wavelength  of  532  nm. 

All  the  rays  entering  the  test  section  are  in  the  z- 
direction.  In  a  flow  with  the  density  varying  as  in  Eq. 
(1),  a  ray  will  be  independent  of  x,  and  will  stay  in 
a  y-z  plane  reducing  the  problem  to  two-dimensions. 
With  reference  to  Fig.  2  for  the  coordinates  (y,z),  a 
ray  in  the  flow  is  then  governed  by  the  equations  that 
are  readily  derived  from  the  eikonal  equation7 


where  s  is  a  one-dimensional  spatial  coordinate  vari¬ 
able  along  ray  path. 


The  ray  direction  is  specified  in  terms  of  the  di¬ 
rection  cosines  dy/ds  and  dz/ds,  which  satisfy 


with  dz/ds  —  0.  Eq.  (2)  is  integrated  as 

dz  _  C 
ds  n  ’ 


(4) 


(5) 


with  C  being  the  integration  constant  to  be  deter¬ 
mined.  All  the  rays  entering  the  test  section  are  in 
the  z-direction,  and  thus 


dz 

ds 


IjsO 


1. 


(6) 


It  follows  then  that  C  =  n(yo),  Vo  being  the  y- 
coordinate  of  the  ray  entering  the  test  section,  and 
then  Eq.  (5)  is  written  as 


dz  _  n(y p) 
ds  n(y)  ’ 

From  Eqs.  (5')  and  (4),  one  obtains 


Division  of  Eq.  (5')  by  Eq.  (7)  yields 
dz  1 


(5') 


(7) 

(8) 


Integration  from  the  ray  entrance  point  (y  =  y0,  z  = 
0)  to  a  point  (y,z)  on  the  ray  path  in  the  test  section 
can  be  written  as 


The  refractive  index  of  air  is  a  function  of  y  only, 
as  given  in  Eq.  (1).  The  air  density  is  0.001226  g/cm3 
under  ambient  conditions,  and  thus  the  term  contain¬ 
ing  the  air  density,  Kp,  is  close  to  0.000278,  and  Eq. 
(9)  can  be,  to  first  order  approximation,  written  as 


^  /S  fly 

ho  \/2 K(p(y)  -p{yo)) 


(10) 


This  equation  will  be  solved  for  a  mathematical  model 
characteristically  representing  a  real  flow  subjected  to 
a  large  density  gradient,  and  solutions  will  be  exam¬ 
ined  to  account  for  interferogram  distortions. 

In  the  interferograms  being  analyzed  in  this 
study,  a  large  interferogram  distortion  appears  pri¬ 
marily  in  the  vicinity  of  the  upper  surface  of  the  air¬ 
foil  near  the  leading  edge  for  conditions  of  high  levels 
of  positive  lift.  In  such  a  flow,  the  air  density  is  the 
smallest  on  the  surface,  increases  with  increasing  dis¬ 
tance  away  from  the  surface,  and  then  asymptotically 
reaches  the  free  stream  density  Poo  .  The  density  gra¬ 
dient  is  positive  and  the  greatest  on  the  surface,  and 
decreases  with  increasing  y.  The  largest  ray  deflec¬ 
tion  occurs  near  the  airfoil  surface  (y=0).  On  the 
other  hand,  ray  deflections  are  expected  to  be  small 
for  a  short  travel  time  in  the  flow.  In  other  words,  the 
y-coordinate  variation  of  a  ray  will  be  small.  Under 
these  conditions,  the  density  can  be  expanded  as 


p{y)  =  P,  +  ay  -  by1  +  •••,  (11) 

where  p,  is  the  density  on  the  airfoil  surface,  and 

dp\ 


a  = 


dy 


(12) 


»=o 
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iff. 

2  dy3 


(13) 


In  such  a  flow  as  described  above,  a  and  b  are  both 
greater  than  0. 


With  Eq.  (11)  inserted  into  Eq.  (10),  the  inte¬ 
gration  yields 

V  =  Pd  +  w  [1  -  cos(/?z)],  for  0  <  z  <  L,  (14) 
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where 


1 


-  =  |  -  y°, 

(15) 

0  =  s/2Kb. 

(16) 

The  tangential  slope  of  the  ray  is  obtained  as  the 
derivative 


Projected  onto  the  image  plane,  this  virtual  displace¬ 
ment  is  translated  as  an  image  displacement,  6,  which 
is  obtained  as 

6(z)  =  m(z)6(z).  (23) 

Here  m  is  the  magnification,  which  depends  on  the 
location  of  the  object  plane  (z). 


j/  =  tv0  sin(/3  z),  for  0  <  z  <  L.  (17) 


Once  a  ray  exits  the  test  section,  it  will  follow 
a  straight  line,  and  is  specified  by  a  linear  equation 
which  is  determined  by  the  y-coordinate,  yi,  and  the 
slope,  of  the  ray  at  exit  (z=L)  of  the  test  section. 
These  parameters  are  obtained  by  inserting  z=L  re¬ 
spectively  into  Eqs.  (14)  and  (17)  as 


Vl  =  Vo  +  w  [1  —  cos(0L)]t 

(18) 

and 

y/l  =  w0sin(PL). 

(19) 

The  equation  for  the  straight  ray  is  then 

V  =  1 4  (*-£)  +  W-- 

(20) 

As  governed  by  this  equation  and  Eq.  (14),  rays  are 
displayed  in  Fig.  2.  As  can  be  seen  readily  from  Eqs. 

ii4)  and  (151,  the  in-flow  ray  deflection  is  the  greatest 
or  the  incident  ray  of  y0  =  0,  and  decreases  with 
increasing  pt>.  From  Eqs.  (15)  and  (17),  the  ray  slope 
at  z=L  is  also  seen  to  be  the  greatest  for  yo  =  0,  and 
to  decrease  with  increasing  yo- 

4.  Image  Displacement 

Eq.  (20)  will  now  be  examined  to  study  imag¬ 
ing.  When  the  image  is  focused  on  an  object  plane 
(z=constant),  an  incident  ray  of  y=  yo  will  constitute 
an  image  of  the  object  that  is  on  the  ray  path  but  will 
appear  virtually  to  be  located  at  (y,z)  with  y  being 
determined  by  Eq.  (20).  Without  the  density  variar 
tion,  the  ray  would  not  deflect  in  the  flow  (as  assumed 
in  deflection-free  interferometry),  and  would  remain 
as  a  straight  line  parallel  to  the  z-axis.  Such  a  ray 
would  be  governed  by  the  equation 

V  =  Vo-  (21) 

Defining  the  virtual  displacement,  6,  as  the  difference 

between  Eqs.  (20)  and  (21),  one  obtains 

$(*)  =  l4  (*  ~  L)  +  VL  -  Vo 
=  w0  sin(P  L)(z  —  £.)  + 

w[\  -  co»(0L)].  (22) 

For  0L  <  1, 

6(z)  9S  (a  -  2byo)KL(z-±).  (22 a) 


We  digress  here  to  briefly  review  the  image  for¬ 
mation  by  a  spherical  mirror,  in  terms  of  three  pa¬ 
rameters:  f,  the  focal  length  of  the  mirror,  and  p  and 
q,  the  distances  from  the  mirror  respectively  to  the 
object  plane  and  the  image  plane.  These  parameters 
are  interdependent  through  the  mirror  formula 


1  1  _  1 
P  +  ?  ~  /’ 
and  the  magnification  is  given  by 


P 


(24) 

(25) 


It  follows  that,  with  an  image  plane  chosen,  the  loca¬ 
tion  of  the  object  plane  is  determined,  and  vice  versa. 


The  virtual  displacements  on  the  object  plane  and 
their  projection  on  the  corresponding  image  planes 
are  shown  in  Fig.  3.  In  the  figure,  the  focal  length 
is  10  in  a  relative  scale  used  for  distances.  The  sub¬ 
scripts  0  and  L  are  used  to  indicate  respectively  the 
object  planes  z=0  and  z=L  When  focused  on  the  en¬ 
trance  plane  (z=0),  po=19  and  q0=21.1;  and  focused 
on  the  exit  plane  (z=L),  p/,=15  and  q/.=30.  For  the 
incident  ray  with  yo=0,  the  virtual  displacements  are 
shown  as  So  and  Si  on  the  respective  object  planes, 
and  their  projections  on  the  image  planes  are  So  and 
Si.  Note  that,  as  evaluated  from  Eq.  (22),  60  <  0 
and  Si  >  0.  As  displayed  in  this  figure,  So  <  0, 
Si  >  0.  It  follows  then  that  the  flow  image  pene¬ 
trates  into  the  image  of  the  airfoil  by  -So  when  focused 
on  the  entrance  plane  z=0,  as  can  be  observed  in  Plate 
III  a  and  b.  On  the  other  hand,  when  focused  on  the 
exit  plane  z=L,  the  image  is  separated  from  the  the 
surface  of  the  airfoil  image  by  SL,  leaving  the  gap  to 
be  completely  dark  as  seen  in  Plate  IV  a  and  b. 

The  analysis  can  be  applied  to  a  case  of  density 
variation  along  the  surface  to  account  for  the  dark 
spot  shown  near  the  leading  edge  of  the  airfoil  as  in 
Plate  II.  In  this  case,  the  y-coordinate  variable  is  re- 

E laced  by  a  coordinate  variable  parallel  to  the  surface 
ut  perpendicular  to  the  s-axis. 


5.  Quantitative  Evaluation  of  Image 
Displacement  for  A  Case  of  Density  Variation 

We  now  consider  a  density  variation  which  is  typi¬ 
cal  of  dynamic  stall,  as  displayed  schematically  in  Fig. 
3.  It  is  mathematically  modeled  by  the  equation 

p(y)  =  p.  +  Ap(l  -*-»'»•).  (26) 

Here  A p  —  Poo  —p$,  and  y,  is  the  y-coordinate  where 
the  air  density  is  dose  to  Poo  with  a  difference  of  A  p/e. 
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This  model  exhibits  the  characteristics  of  a  flow  in  the 
vicinity  of  a  surface  such  as  the  upper  surface  of  the 
airfoil  near  the  leading  edge  for  cases  of  high  levels  of 
positive  lift.  The  air  density  increases  with  increasing 
y  and  asymptotically  reaches  the  free  stream  density 
poo .  The  density  gradient  is  the  greatest  on  the  sur¬ 
face  and  decreases  with  increasing  y.  For  this  flow 
model,  one  readily  obtains,  on  inserting  Eq.  (26)  into 
Eqs.  (12)  and  (13), 


- 

Vt  ' 


(12a) 


iA/> 


(13a) 


On  inserting  these  into  Eqs.  (15)  and  (16),  one  ob¬ 
tains 


v>  =  Vt  -  Vo, 
-  =  0-0167^ 
V* 


(15a) 

(16a) 


where  y  =  Ap/po ,  and  Kpo  =  0.000278  has  been 
used. 


As  an  example,  we  consider  an  interferogram 
taken  with  M=0.3  and  a=16°.  By  using  the  fringe 
counting  method,  we  determined  7  9J  0.30,  and 
y«  SE  0.52cm.  For  yo  =  0,  w  =  0.52cm  and  0 
=  0.0176  cm-1.  Using  Eq.  (22),  one  obtains  bL  — 
0.05  cm.  In  this  case,  m=2.77;  thus  the  image  dis¬ 
placement  61  =  0.14  cm.  Focused  on  the  plane  z  =  0, 
60  =  -  0.049,  m=2.22,  and  6 0  =  -  0.11  cm.  The  inter¬ 
ferogram  was  taken  with  the  system  focused  on  three 
triangular  objects  patched  on  the  entrance  window 
(z=0j.  The  fringe  penetration  into  the  airfoil  image 
was  measured  from  a  direct  comparison  of  the  inter¬ 
ferogram  with  the  airfoil  image  taken  without  flow. 
The  measurement  was  0.112  cm,  and  is  in  an  excel¬ 
lent  agreement  with  the  calculated  value  of  -Sq,  0.11 
cm. 


6.  Intensity  of  Deflected  Light 

The  rays  are  conserved,  and  since  the  rays  are  al¬ 
most  parallel  to  the  z-axis,  the  ray  density  is  inversely 
proportional  to  the  vertical  separation  dy  of  two  ad¬ 
jacent  rays.  The  light  intensity  which  is  proportional 
to  the  ray  density  is  then  subjected  to  the  relation 
I(y,z)dy  =  Iodyo,  and  thus 

/(»,*)  =  Jj;,  (27) 


Here  I(y,z)  is  the  light  intensity  at  (y,x),  and  Io  is  the 
incident  light  intensity  at  *=0,  which  is  assumed  to 
be  uniform.  From  Eqs.  (27)  and  (14),  one  obtains 


I  =  — — — . 
eo$(Pz) 


(28) 


This  equation  is  valid  only  in  the  region  specified  by 

^[1  -cos(0z)]  <  y  <  a/26.  (29) 

Below  the  lower  limit,  there  is  no  light,  and  the  region 
is  completely  dark.  As  y  approaches  the  upper  limit, 
the  expansion  in  Eq.  (11)  loses  its  validity.  Within 
the  valid  region,  Eq.  (28)  shows  that  the  light  in¬ 
tensity  increases  slightly  with  increasing  z,  but  does 
not  depend  on  y.  The  present  analysis  generates  a 
shadowgraph  as  a  dark  spot;  however,  the  intensity 
of  the  displaced  image  hardly  fluctuates.  This  is  ex¬ 
pected  because  the  second  derivative  of  the  air  density 
is  constant  as  seen  from  Eq.  (11).  The  present  analy¬ 
sis  produces  also  schlieren  effects  as  expected  from  Eq. 
(19).  Although  the  deflection  angle  depends  implic¬ 
itly  on  y  via  yo ,  the  intensity  of  the  schlieren  images 
is  again  constant.  The  density  variation  around  the 
airfoil  is  a  well-behaved  field;  thus  intensity  fluctua¬ 
tion  of  shadowgraphs  or  of  schlieren  images  will  be 
minimal. 


7.  Concluding  Remarks 

Detailed  understanding  and  correction  of  inter¬ 
ferogram  distortion  is  critical  to  accurate  analysis  of 
the  compressible  flow  that  occurs  on  airfoils  at  high 
angles  of  attack.  Correction  of  these  distortions  is 
absolutely  necessary  if  any  form  of  flow  control  is  to 
be  considered  in  tins  unsteady  compressible  flow  en¬ 
vironment. 

An  analytical  method  was  formulated  to  account 
for  distortions  of  interferograms  used  in  studies  of  dy¬ 
namic  stall  flow.  The  analysis  correctly  predicts  the 
distortions.  However,  since  the  method  assumes  a  one 
dimensional  density  variation,  it  should  be  expanded 
and  combined  probably  with  a  numerical  scheme  to 
completely  correct  the  interferogram  distortions. 
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Fig.  3  Virtual  and  image  displacements  due  to  ray  bending  in  flow. 
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Fig.  4  Air  density  vs.  y  for  a  model  case. 
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Plate  I  Interferogram  distortion  near  leading  edge. 


Plate  II  Interferogram  distortion:  A  dark  spot  near  leading  edge. 
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Plate  III  Image  focused  on  entrance  plane:  a=10#;  a)  M=0.4,  b)  M=0. 
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